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ABSTRACT:  Ab  initio  calculations  have,  been  performed  to  determine  potential 

2  2  + 

curves  and  molecular  properties  for  the  X  n  and  A  Z  states  of  the  LiO  molecule. 

The  calculations  use  a  conventional  configuration  interaction  (Cl)  method  in 

which  the  Hartree-Fock  configuration  is  taken  as  reference  configuration  and 

only  valence  shells  are  correlated.  The  molecular  orbital  (MO)  basis  set 

used  in  the  Cl  calculations  is  composed  of  the  Hartree-Fock  orbitals  and 

additional  MO's.  These  additional  MO's  are  formed  by  truncating  a  set  of 

pseudonatural  orbitals  (PSNO's)  obtained  as  the  natural  orbitals  of  a  Cl 

calculation  on  a  single  pair  of  valence  electrons.  The  main  results  are 

R  =  1,695^,  D°  =  3.37  eV,w  -  851.5  cm”1  (7Li160) ,  y  -  6-76  Debye  for  the 
e  e  e 

X2n  state  and  R  =  1.599A,  D°  «  4.90  eV,w  =  866.8  cm”1  (?Li  160;  ,  y  =  5.96 

6  6  C 
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Debye  Cor  the  A  Z  state.  The  computed  Dq(X  H)  ■  3.37  eV  is  in  good 

agreement  with  the  observed  value  of  3.39  ±  0.26  eV.  The  other  results 

are  also  believed  to  be  accurate  to  within  a  few  percent.  The  computed  term 
2  •>  +  -1 

energy  Tfi(X  IT  *►  A*T  )  is  2330  cm  with  the  Cl  wave  functions,  in  surprising 
agreement  with  the  Hartre^-Fock  value  of  2342  cm 
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ABSTRACT:  Ab  initio  calculations  have  been  performed  to  determine  potential 
curves  and  molecular  properties  for  the  X2H  and  A25f  states  of  the  LiO  molecule. 
The  calculations  use  a  conventional  configuration  interaction  (Cl)  method  in 
which  the  Hartree-Fock  configuration  is  taken  as  reference  configuration  and 
only  valence  shells  are  correlated.  The  molecular  orbital  (MO)  basis  set 
used  in  the  Cl  calculations  is  composed  of  the  Hartree-Fock  orbitals  and 
additional  MO's.  These  additional  MO's  are  foi'med  by  truncating  a  set  of 
pseudonatural  orbitals  (PSNO's)  obtained  as  the  natural  orbitals  of  a  Cl 
calculation  on  a  single  pair  of  valence  electrons.  The  main  results  are 

R  =  1.695A,  D°  =  3.37  eV,we  =  851.5  cm"1  (?Li160),  y  =  6.76  Debye  for  the 
X^n  state  and  Rp  =  1.599A,  D°  =  A. 90  eV,u  =  866.8  cm"1  (7Li160),  y  =  5.96 

5  +  e  e  o 

Debye  for  the  A^E  state.  The  computed  D°(X^II)  =  3.37  eV  is  in  good  agreement 


with  the  abserved  value  of  3.39  ±  0.26  eV.  The  other  results  are  also 
believed  to  be  accurate  to  within  a  few  percent.  The  computed  term  energy 
Te  (X2n  -+  A2Z+)  is  2330  cm”-*-  with  the  Cl  wave  functions,  in  surprising  agree¬ 
ment  with  the  Hartree-Fock  value  of  2342  cm-^. 
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INTRODUCTION 

Gaseous  nonoxides  are  important  high  temperature  systems  being  the 

subject  of  numerous  experimental  investigations  in  recent  years. *  From  a 

theoretical  point  of  view,  the  LiO  molecule  is  one  of  the  smallest  and 

simplest  gaseous  monoxides;  however,  very  few  experimental  data  are  presently 

available,  in  fact,  the  dissociation  energy  is  only  property  for  which  a 

reliable  value  is  known  experimentally.  There  is  a  measurement  of  the 

vibrational  frequency  v(l-O)  of  LiO  in  an  inert  gas  matrix  but  its  gas 

2 

phase  value  is  not  known.  It  has  been  established  that  there  is  a 

considerable  shift  in  the  value  of  v  measured  in  an  inert  gas  matrix 

from  its  gas  phase  value  for  molecules  with  small  mass  and  large  dipole 

3 

moment,  in  particular  for  the  liF  molecule.  Since  the  LiO  molecule  is 

similar  to  the  LiF  molecule,  one  would  also  expect  a  considerable 

difference  between  v  measured  in  an  inert  gas  matrix  and  in  the  gas 

« 

phase. 

In  this  work,  we  have  determined  potential  curves  and  other  properties 

2  2  + 

for  the  two  lowest  states,  X  n  and  A  £  i  of  the  LiO  molecule.  These  results 

have  been  obtained  with  large  Cl  wavefunctions  including  ^1500  symmetrized 

configurations.  Techniques  of  quantum  mechanical  calculations  have  been 

improved  recently  and  it  is  now  possible  to  calculate  with  reasonable 

accuracy  many  properties  of  system  like  LiO.  In  other  words,  the  ab  initio 

results  may  be  used  to  fill  gaps  in  experimental  data  and  that  is  the  main 
"  .  *'  ~  Continue 

purpose  of  this  work.  c  ,, 

r  r  rutlOW  l  i  j  .  '</,  ■  '  I 

■  It  is  still  impracticable  to  carry  out  Cl  calculations  which  give 

potential  curves  with  small  absolute  errors  even  for  a  system  of  the  size 
of  LiO.  However,  it  is  only  necessary  to  obtain  a  potential  curve  which 
is  closely  parallel  to  the  actual  curve  in  order  to  derive  reliable 

■  / /; ’/)’//  iJll  I'J'tt  ,\C  't\! 
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spectroscopic  constants  and  other  properties.  In  the  present  calculations, 
we  have  only  considered  correlation  effects  due  to  interactions  of  the 
valence  shell  electrons  with  themselves.  If  we  start  from  llartree-Fock  or 
near  Hartrce-Fock  solutions,  we  expect  that  the  contribution  to  the  total 
correlation  energy  of  interactions  of  the  core  electrons  (Is  in  Li  and  0) 
with  themselves  and  with  the  valence  shell  electrons  will  remain  constant 
along  the  potential  curve.  We  also  expect  that  the  neglect  of  these 
interactions  will  have  a  small  effect  on  other  properties.  In  practice, 
this  restriction  allows  us  to  treat  the  valence  shell  correlation  more 
accurately  and  in  a  more  uniform  manner  along  the  potential  curve.  This 
is  possible  because  we  are  able  to  include  more  valence  shell  correlating 
orbitals  and  configurations  in  the  Cl  calculation  than  would  be  possible 
if  all  electrons  were  correlated. 

In  Cl  calculations  it  is  of  utmost  importance  to  have  a  proper  orbital 
basis  and  proper  choice  of  configurations.  To  meet  this  end,  the  following 
procedure  was  used  for  the  present  calculations:  1)  Obtain  Hartree-Fock 
(HF)  wavefunctions  for  various  intemuclear  distances;  2)  Obtain 

4 

pseudonatural  orbitals  (PSNO's)  for  an  appropriate  pair  of  electrons; 

3)  Form  the  MO  basis  for  the  Cl  calculation  by  adding  a  truncated  set  of 
the  PSNO's  to  the  Hartree-Fock  occupied  orbitals;  4)  Perform  the  Cl 
calculation  using  the  Hartree-Fock  configuration  as  the  reference  configura¬ 
tion  and  including  all  configurations  generated  by  single  and  double 
substitution  of  the  valence  orbitals  into  all  the  PSNO's  used.  The  detail 
of  the  procedure  will  be  given  in  the  next  section. 

This  procedure  is  necessarily  pragmatic  in  its  nature.  In  addition 
to  assessing  the  quality  of  the  computed  wavefunctions  and  properties, 
comparison  of  the  results  with  known  experimental  data  serves  the  purpose 

..  L  i  : 
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of  testing  the  general  applicability  of  the  procedure.  This  is  another 
important  purpose  of  this  work. 

These  calculations  arc  performed  with  the  ALCHEMY  system  of. programs 
developed  by  the  quantum  chemistry  group,  IBM  Research  Laboratory,  San  Jose. 


CALCULATIONS 

A  basis  set  of  Slater-type  functions  is  constructed  based  on  the 
accurate  atomic  basis  sets  of  P.  S.  Bagus  et  al. for  both  Li  and  0.  The 
atomic  basis  sets  are  supplemented  with  polarization  functions  needed  for 
the  HF  calculation  and  also  other  Slater-type  functions  of  higher  angular 
momentum  needed  for  describing  axial  correlation  effects.  The  number  and 
exponents  of  these  supplementary  functions  are  chosen  so  that  the  basis 
set  will  span  an  adequate  space  without  performing  optimization  of 
exponents.  The  result  is  a  rather  large  basis  set  of  36a,  22tt,  106,  4$ 
type  functions,  which  is  listed  in  Table  I. 

2  2  2  2  3 

With  this  set,  single  configuration  SCF  calculations,  lo  2o  3o  4a  lit 

2  2  2  2  4  2  +  . 

for  the  X  n  state  and  la  2a  3a  4o1ti  for  the  A  E  state,  are  performed  at  a 

number  of  points  around  the  equilibrium  intemuclear  distance  R^.  The 

same  basis  set  is  used  for  both  states  and  for  all  separations.  These 

single  configurations  do  not  dissociate  to  the  proper  states  of  the 

separated  atoms  for  either  the  X  H  or  A  I  state.  However,  for  the  range 

of  R  (around  Re)  of  interest,  we  believe  that  the  single  Hartree-Fock 

r  ■  *  I  I* 

configuration  is  sufficiently  dominant  to  be  suitable  as  the  "reference 
state"  in  subsequent  Cl  calculations.  The  configurations  in  the  Cl 
calculations  are  constructed  and  described  as  zero,  one,  two,  etc. 
replacements  from  the  reference  configuration(s) .  If  several  configurations 

7 . 
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huvo  largo  weights  In  a  Cl  waveiunctiuA,  it  nay  bo  important  to  inrludo 

tlicoo  configurations  as  part  of  the  "reference  state."  We  have  verified 

that  this  is  not  the  case  in  our  calculations  by  examining  the  weights  of 

configurations  and  the  occupation  nuabers  of  the  natural  orbitals. 

The  results  of  SCF  calculations  indicate  that  our  choice  of  the  basis 

set  is  a  reasonable  one  so  far  as  the  total  SCF  energy  is  concerned. 

Further,  there  were  no  problems  of  basis  set  linear  dependence  for  the  SCF 

calculations.  However,  ti.s  basis  set  was  redundant  for  the  Cl  calculations. 

Tnis  redundancy  is  manifested  by  several  virtual  Hartree-Fock  orbitals  with 

large  Cv50-100)  positive  and  negative  coefficients,  causing  a  substantial 

loss  of  accuracy  in  the  four  index  transformation  of  electron  repulsion 

Integrals  from  integrals  over  STO's  to  integrals  over  molecular  orbitals. 

It  was  thus  necessary  to  reduce  the  space  spanned  by  the  orthogonal 

complement  of  the  occupied  Hartree-Fock  orbitals.  This  was  done  in  the 

% 

following  way.  First,  the  overlap  matrix  of  tie  STO  basis  functions  was 

diagonalized.  The  eigenvectors  associated  with  eigenvalues  less  than  a 

certain  threshold  were  deleted.  The  SCF  calculation  was  then  performed 

in  the  space  spanned  by  the  remaining  eigenvectors  of  the  overlap  matrix. 

(The  SCF  energies  for  LiO  of  these  calculations  differed  only  in  the  fifth 

decimal  place  from  those  obtained  using  the  full  STO  basis  set.)  The 

virtual  Hartree-Fock  orbitals  of  this  new  calculation  are  well  behaved,^ 

and  were  orthogonalized  to  the  occupied  Hartree-Fock  orbitals  obtained 

without  truncating  the  space  spanned  by  the  STO  basis  set.  (This 

orthogonalization  introduced  only  very  small  changes  in  the  virtual 

orbitals.)  These  orbitals  plus  a  set  of  orthogonalized  6  and  $  type  STO's 
♦ 

formed  a  molecular  orbital  basis  set  for  PSNO  calculations;  the  sets  are 
2  2  + 

denoted  by  II)  and  T.  ).  The  sets  <j>^  consisted  of  31o ,  20ir,  106, 

i,..1 7  /  //,/■>  [j •  •: » [*7 ?  iV'  i  f  '  •’  .7 
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and  4:  typo  orbitals  for  2.3  0.1  a.u.  and  32o,  20*,  10*,  and  4*  tvpa 

for  3.2  <_  R  <  3.9  a.u. 

Since  L10  ie  Ionic  in  the  vicinity  of  ita  equillbriun  lntemuclear 
separation  (R  )  and  the  correlation  effect  of  valence  shell  electrons  alone 
is  consldeied  here,  the  problem  reduces  essentially  to  correlating  the 
valence  electrons  in  0**.  Weiss'  work2  on  the  first  row  atowi  shows  that 
the  psoudonatural  orbital  irons forest  ion  on  any  single  pair  of  valence 
electrons  yields  a  rapidly  converging  and  varlatlonally  optima  set  of 
virtual  orbitals.  We  have  calculated  several  of  the  pair  correlation 

i  2  + 

energies  for  the  valence  shells  of  both  the  * H  and  A  £  states;  angular 

Halts  as  well  as  total  pair  correlation  energies  have  been  obtained.  Wr 

hav'  chosen  the  1 — 1'  pair  for  the  calculation  of  all  but  one  set  of 

PSHO's.  This  was  done  because  the  pair  correlation  energy,  sore  precisely 

the  appropriate  angular  Halt  of  the  pair  correlation  energy,  vas  largest 

for  this  pair.  The  one  exception  is  the  set  of  o  PSNO's  for  the  j^n  state; 

here  the  4o-4o  pair  has  been  used  since  the  e-llalt  of  the  pair  correlation 

energy  for  this  pair  is  larger  than  for  the  ls-lrr  pa!**. 

The  configurations  along  vlth  the  coupling  schemes  used  in  the  PSNO 

calculations  arc  listed  below: 

2 

For  the  II  state 


(lo22o23o2lr3)2n(iMno)1£+  , 

j  /. *  *7  i  i  *  •  (/  /  />•  i  iff 

j  2  2  2  2  2  1  +  Fn  :.  t  t :  ' 

<1cj  2<T3cj  4c>  litrn(mXnXrr  ,  X  4  o  , 


2  + 

and  for  the  l  state 


•.  P  Lin  v  •* 


2»  2,,  2.  ,  2.2..+,  .  .  .1,.+ 
(lo  2o  3a  4a In  )  1  (mXnX)  , 

it  .  .  ,/./•  />; "/.  J.  !  * 
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where  m,n  ■  4,. ..,31  for  >  *  0,  w»n  *  2,.. .,20  for  *  «  1, 

n,n  ■  1 , ....  10  for  *  **  2,  and  n,n  »  1,...,4,  for  >  •  3. 

Out  of  the  full  set  of  PSSO's,  llo,  10",  55  and  2:  virtual  orbitals, 
all  of  which  have  the  occupation  number  greater  than  )0~^,  arc  selected 
and  corJbined  with  the  Hartrce-)‘ock  orbitals  to  produce  orthonorralized 

2  2  * 

orbital  seta,  ;^(  F.)  and  {g(  1).  The  composition  of  i-g  is  15*?,  llr, 

5&  and  2c  type  orbitals.  The  effect  on  the  total  energy  of  truncation 
in  this  net hod  of  selecting  virtual  orbitals  has  been  studied  and  will  be 
discussed  later  in  this  section. 

Configuration  wavefunctlons  are  constructed  with  the  orbital  set 
Cg.  All  configurations  are  Included  which  have  formally  non-vanishing 
Hamiltonian  matrix  elements  with  the  single  HF  reference  configuration, 
and  which  could  be  constructed  by  single  and  double  excitations  from  the 
valence  shell  orbitals  in  the  reference  configuration  (in  this  system, 
the  3o,  4o,  and  la  orbitals)  into  all  virtual  orbitals.  We  have  Included 
all  single  excitations  even  though  many  of  them  have  numerically  zero 
matrix  elements  with  the  reference  configuration  because  of  Brillouin's 

g 

theorem.  This  was  done  to  ensure  a  proper  wavofunction  for  the  evaluation 

9 

of  one  electron  property  such  as  the  dipole  moment. 

It  should  be  noted  thet  for  some  orbital  configurations  which 

result  from  double  excitations,  the  orbital  and  spin  angular  momentum 

coupling  give  more  than  two  linearly  independent  states  of  the  dbslrcd 

2  2  2 

total  symmetry.  [For  example,  a  Is  -  na  excitation  for  the  n  state 
gives  rise  to  three  states:  ' 

(1)  4o2l*(2n)nr2(V)  , 

(2)  4n2l«(2U)nn2(1&)  , 

(3)  4o2M2!!)nvVi:+)  .) 

.  .  i  i.>*  i 
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All  of  these  states  nay  have  non -vanish tup  llAniltonlon  matrix  elements 

with  the  reference  configuration.  However*  these  matrix  elenentn  are  all 

expressed  in  ternt  of  only  two  different  two-electron  integrals.,  [For 

o  2 

example*  above  they  are  K  (ls*li)  and  K  (ln*lr).]  For  these  cases  a 

reduction  of  the  number  of  configurations  can  be  accomplished  by  taking 

proper  linear  combinations  of  these  configurations  so  that  only  two 

resultant  configurations  have  nor.- vanishing  Hamiltonian  matrix  elements 

with  the  reference  configuration.  This  is  dona  whenever  applicable  in  the 

present  calculations  resulting  in  a  considerable  reduction  of  the  number 

of  configurations  to  be  included.  The  actual  numbers  of  configurations  used 

2  2  ♦ 

arc  1632  and  1A58  for  the  X  I!  and  A  Z  states,  respectively.  The  some  set 

of  configurations  is  used  for  all  values  of  R  except  R  *  SO  a.u. 

In  order  to  estimate  the  dissociation  energy,  a  similar  calculation 
*> 

Cor  the  X"ll  state,  but  with  a  different  reference  configuration,  was 
carried  out  at  R  •  SO  a.u.  where  for  practical  purposes  the  system  is  in 

the  separated  atom  limit.  The  reference  configuration  used  here  is 

2  2  2  3  2  3 

lo  2o  3a  A~5cl"  which  dissociates  properly  to  Li:  S  and  0:  P.  At  this 

distance  all  orbitals  are  necessarily  atomic  like  and  atomic  symmetry  can 

be  taken  into  account;  this  reduces  the  number  of  configurations  for  a 

given  set  of  orbitals.  To  be  consistcr.*:  with  the  calculations  at  other 

values  of  R  it  is  only  necessary  to  consider  the  correlation  effect  of 

the  valence  electrons  in  the  oxygen  atom.  We  used  the  llartree-Fock 

orbitals  plus  all  virtual  orbitals  belonging  to  the  oxygen  atom;  thus  the 
2 

set  <?^(  II)  at  R  -  50  a.u.  consists  of  18a,  10tt,  56  and  2$  type  orbitals. 

2 

A  PSN'O  calculation  was  not  carried  out  here  since  the  set  4>^(  n)  is  small 
enough  that  all  orbitals  cun  be  included  in  the  final  Cl.  The  error  due 
to  the  truncation  of  orbitals  another  values  of  R  turned  out  to  be 

...»  t )  '  it  .  *  . 


8 


negligible  for  estimating  the  dissociation  energy  (see  the  discussion 
below).  The  configurations  for  the  final  Cl  calculation  at  R  =  50  a.u. 
were  selected  in  a  manner  similar  to  that  described  above  except  that  the 
set  4^  is  used  instead  of  4>g.  However,  for  the  cases  where  more  than 
two  states  can  be  obtained  by  angular  momentum  coupling  from  a  given 
orbital  configuration,  the  reduction  of  the  number  of  these  states  was 
not  completely  carried  out-  The  total  number  of  configurations  used  was 
1431. 

As  mentioned  above,  the  effect  of  the  truncation  of  orbitals  was 

2 

studied  at  R  =  3.2  a.u.  for  the  n  state  by  computing  pair  energies 

2  2 

separately  with  orbital  sets  <J>.(  II)  and  $_(  11).  The  pair  energy  is 

A  D 

defined  here  as  the  energy  improvement  obtained  in  a  Cl  calculation  using 

single  and  double  excitation  from  one  pair  of  HF  space  orbitals  at  a  time. 

The  results  are  shown  in  Table  II.  Truncation  errors  are  given  in 

magnitude  and  in  percentage.  As  expected,  the  3o  pair  energy  is  the  worst 

case  and  gives  a  3%  error.  The  error  in  the  sura  of  all  pair  energies  is, 

however,  1.4%  indicating  our  method  of  selecting  orbitals  is  satisfactory. 

2  + 

No  investigation  was  made  for  the  E  state  but  it  is  expected  that  the 

2 

magnitude  of  the  errors  should  be  the  same  as  for  the  n  state. 


RESULTS  AND  DISCUSSION 


2  2  + 

The  SCF  and  Cl  results  for  the  X  n  and  A  E  states  of  LiO  computed  at 
intemuclear  distances  very  close  to  their  respective  equilibrium  inter- 
nuclear  distances  are  listed  in  Table  III.  The  Cl  energies,  the  energy 
improvements  in  Cl,  and  selected  molecular  properties  computed  at 'various 
intemuclear  distances  are  shown  in  Tables  IV  and  V  for  the  X^n  and  A  E 


states,  respectively. 


The  Cl  potential 


tentii 


■  U.  j’iii  ■  l  JO 


curves  are  also  shown  in 


Figure  1. 
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The  previously  calculated  SCF  energy  at  R  =  3.184  a.u.  for  the  X“H  state 
is  -82.3111  a.u.^  and  our  value  is  -82.3115  a.u.,  an  improvement  of 
.0004  a.u.  The  SCF  energies  obtained  in  the  present  calculations  should 
be  close  to  the  exact  Hartree-Fock  values;  at  least  within  .0001  a.u. 

Our  SCF  wavefunctions  satisfy  the  virial  theorem  very  well  and  the  computed 
net  forces  on  the  nuclei  are  very  small.  These  observations  also  support 
our  claim  of  having  reached  the  Hartree-Fock  limit.  The  same  observations 
also  hold  for  the  Cl  wavefunctions  although  the  actual  agreements  are 
slightly  worse  than  those  for  the  SCF  wavefunctions. 

2 

The  non-relativistic  total  energy  of  LiO  in  the  X  )1  state  is  about 
-82.666  a.u.,  estimated  from  the  experimental  dissociation  energy  of 
.1247  a.u.l  and  atomic  energy  data.^  The  correlation  energy  of  LiO  is 
-.355  a.u.  out  of  which  we  estimated  that  -.250  a.u.  is  the  valence  electron 

correlation  energy  (EyC) •  This  estimate  for  E^c  is  obtained  as  follows: 

12 

first  taking  Eyc  of  atomic  oxygen  to  be  -.193  a.u.,  the  rest  of  the 
correlation  energy  for  the  separated  atom  system  is  determined  to  be 
-.115  a.u.,  and  assuming  that  this  part  of  the  correlation  energy  remains 
the  same  in  the  molecular  system,  of  LiO  is  obtained.  The  energy 

improvement  obtained  at  is  -.223  a.u.  for  both  states.  This  is  89% 

of  Eyc  and  63%  of  the  total  correlation  energy.  It  may  be  of  interest 
to  point  out  that  the  sum  of  pair  energies  shown  in  Table  II  overestimates 
EyC  by  about  5%.  The  molecular  properties  presented  in  Table  III  show 
remarkable  agreement  between  the  SCF  and  Cl  results.  For  example,  the 
change  in  dipole  moment  is  less  than  2%  indicating  that  there  is  practically 
no  shift  of  the  charges  from  the  SCF  to  the  Cl  wavefunctions.  The  changes 
in  the  expectation  values  of  electronic  coordinates  ore  also  small;  the 

,  -15- 

largest  is  the  1%  change  in  \  p  }  .  This  is  also  reflected  in  a  relatively 
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small  change  of  .059  a.u.  in  the  nuclear  attraction  energy.  Since  the 
change  in  the  potential  energy  is  -.426  a.u.  and  the  kinetic  energy  is 
.203  a.u.,  the  energy  improvement  in  Cl  come  essentially  from  a. sub¬ 
stantial  decrease  in  1/r^  terms. 

2 

The  X  II  state  of  LiO  dissociates  to  the  ground  states  of  separated 
2  3  2  + 

atoms,  namely  to  Li :  S  +  0:  P;  the  A  E  state  dissociates  to  excited  state 

2  3  2  + 

Li:  P  and  ground  state  0:  P.  Another  possible  dissociation  for  the  A  E 

2  1 

state  is  to  Li:  S  +  0:  D;  however  the  atomic  transition  energy  for 

2  2  13  3  1 

Li:  S  -►  Li:  P,  .0679  a.u.,  is  less  than  that  for  0:  P  -+  0:  D, 

13  2 

.0719  a.u.  For  the  X  II  state  SCF  and  Cl  calculations  have  been  performed 

2  + 

at  R  e  50  a.u.  For  the  A  E  state  the  energy  at  R  -  50  a.u.  was  obtained 

14 

by  adding  the  Hartree-Fock  transition  energy,  .0677  a.u.  ,  for  the 

2  2  2 
S  ■+  P  transition  of  atomic  Li  to  the  SCF  and  Cl  energies  of  the  X  n  state 

at  this  same  separation  of  50  a.u.  The  results .are  shown  in  Table  VI 

together  with  computed  dissociation  energies  and  transition  energies  for 

2  2+  2 
X  II  +  A  E.  The  dissociation  energy  De  for  the  X  II  state  computed  using 

Cl  wave functions ,  3.37  eV,  is  in  excellent  agreement  with  the  experimental 

value  of  3.39  i.  .26  eV.^  This  almost  exact  agreement  may  be  fortuitous 

since  we  did  not  get  all  of  the  correlation  energy.  However,  the  good 

2  + 

agreement  strongly  suggests  (1)  that  the  computed  D£  for  the  A  E  state 
should  also  be  quite  accurate,  and  (2)  that  the  procedure  used  in  the 
present  calculations  is  successful  in  giving  wavefunctions  of  the  same 
level  of  accuracy  along  the  potential  curve.  The  analysis  of  the  computed 
potential  curves,  which  will  be  given  below,  also  supports  the  latter 

4 

assertion. 

Table  VII  gives  spectroscopic  constants  calculated  by  a  Dunham 

•  N~ 

analysis  of  the  Cl  and  HF  curves.  The  available  experimental  data  are  the 
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vibrational  frequency  v(l-O)  =  745  cm  *  for  \i^0  and  789  cm  *  for 

6  16  2  2 
LiAD0  measured  in  a  krypton  matrix  for  the  X  n  state.  It  is  known  that 

there  is  a  shift  in  v  measured  in  an  inert  gas  matrix  from  its  gas  phase 

value  particularly  for  molecules  with  small  mass  and  large  dipale  moment. 

For  example,  for  LiF,  which  is  similar  to  LiO  in  mass  and  dipole  moment, 

the  value  of  v  measured  in  a  krypton  matrix  is  820  cm  *  while  its  gas 

phase  value  is  890  cm  \  with  a  shift  of  70  cm  The  computed  values 

2 

of  v(l-0)  for  the  X  II  state  of  LiO,  using  the  Cl  wave  functions ,  are 

825  cm  *  and  870  cm  ^  for  \i  and  ^Li,  respectively.  These  values  are 

shifted  from  those  measured  in  a  krypton  matrix  by  80  cm  *  and  81  cm 

respectively.  Considering  that  LiO  has  less  mass  and  a  slightly  larger 

dipole  moment  (6.76  Debye)  than  LiF  (6.33  Debye), ^  this  shift  of  80  cm  ^ 

is  reasonable.  Thus,  at  least,  the  computed  value  of  u>e  should  be 

accurate  to  within  a  few  wave  numbers.  The  other  higher  derivatives  are 

inevitably  less  accurate  but  they  also  should  not  be  too  far  from  the  true 

values.  It  is  of  interest  to  note  that  the  HF  results  for  spectroscopic 

constants  shown  in  Table  VII  are  in  reasonable  agreement  with  the  Cl 

results  which  are  believed  to  be  very  close  to  the  true  values.  This 

agreement  is  rather  surprising  since  the  HF  results  for  other  ionic 

molecules,  LiF  and  BeO  for  example,  are  usually  significantly  poorer.^ 

The  Cl  value  of  Rfi  is  slightly  larger  than  the  SCF  value.  Our  Cl  value 

2 

of  1.695A  for  the  X  n  state  differs  considerably  from  the  value  of 

Rg  »  1.62A  which  has  been  estimated  from  experimental  data.  The  following 

observations  can  be  made  in  assessing  the  accuracy  of  computed  . 

Firstly,  the  value  of  Re  =  1.62 A  may  not  be  reliable  since  it  is* estimated 

by  an  empirical  formula  using  a  value  of  the  force  constant  obtained  by  a 

2  '  / 1> 

matrix  isolation  technique.  The  second  observation  is  that  for  LiF  and 


BeO  molecules  the  HF  values  of  R  are  1.555A  and  1.29A  and  the  observed 

e 

values  are  1.564A  and  1.33A,  respectively.^  Since  LiO  is  similar  to 

these  molecules,  its  HF  value  is  expected  to  be  smaller  than  the,  actual 

value.  Thus,  the  Cl  value  of  R  should  be  reliable  within  at  least 

e 

.OlA. 

As  mentioned  in  a  previous  section,  because  of  the  method  employed 
for  deleting  unwanted  virtual  orbitals,  the  number  of  o  type  virtual 
orbitals  used  in  the  PSNO  calculations  is  different  for  R  <_  3.2  and 
R  >  3.2  a. u.  The  effect  of  this  difference  on  the  computed  results  is 

generally  negligibly  small.  However,  the  net  force  on  the  nuclei, 

+  51.  listed  in  Tables  IV  and  V  exhibits  a  small  but  distinct 

JjI  0 

discontinuity  at  R  =  3.2  a.u.  Therefore,  other  means  of  deleting  unwanted 
virtual  orbitals  should  be  employed.  Perhaps  an  even  better  way  would  be 
to  avoid  deletion  of  orbitals  altogether  by  carefully  selecting  basis 
functions  with  less  linear  dependence  among  them. 


CONCLUSIONS 

The  procedure  used  in  the  present  calculations  appears  to  work  well 
in  predicting  accurate  values  for  dissociation  energies,  a  property  which 
is  very  sensitive  to  the  correlation  correction.  For  other  properties, 
the  correlation  correction  appears  to  be  small  and  the  results  obtained 

are  believed  to  be  quite  accurate.  Since  the  procedure  can  be  applied 

~  . .  ~  ~  •  (  > 

to  larger  and  more  complex  systems  at  roughly  the  same  level  of  the 
computational  effort  as  that  required  for  LiO,  it  will  be  worthwhile  to 
investigate  further  the  applicability  of  the  procedure  to  more  complex 
systems  than  the  one  studied  in  the  present  work. 
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Table  1.  Slater-type  functions  used  in  SCF  and  Cl  calculations 


Exponent 

n£m 

on  Li 

on  0 

Iso 

4.699,  2.478 

12.418,  6.995 

2  so 

1.500,  .810,  .643 

3.900,  2.922,  1.818 

3so 

2.681,  .500 

8.681,  1.600 

2po 

3.900,  3.000,  2.111, 

8.450,  4.700,  3.744, 

1.076,  .736 

2.121,  1.318 

3po 

.600 

1.300 

3do 

3.000,  2.000,  1.000 

4.000,  3.000,  2.000 

4  fa 

4.000,  2.000 

4.000,  2.000 

2pir 

3.900,  3.000,  2.111, 

8.450,  4.700,  3.744, 

1.076,  .736 

2.121,  1.318 

3prr 

.600 

1.300 

3dTT 

3.000,  2.000,  1.000 

4.000,  3.000,  2.000 

4  fir 

4.000,  2.000 

4.000,  2.000 

3d6 

3.500,  2.500,  1.500 

4.500,  3.500,  2.500 

4f6 

4.000,  2.000 

4.000,  2.000 

4f$ 

4.000,  2.000 

4.500,  2.500 

1  This  Li 
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2 

Table  II.  Comparison  of  pair  energies  for  the  n  state  of  LiO 
computed  with  different  size  orbital  sets  at 
R  =  3.2  a.u.  Energies  are  in  a.u. 


Orbital 

Set 

Difference 

Pair 

A 

B 

Magnitude 

Percentage 

3o3a 

.01587 

.01540 

.00047 

3.0 

3o4o 

.03655 

.03580 

.00075 

2.0 

4o4o 

.02740 

.02735 

.00005 

•2 

3o1tt 

.05723 

.05570 

.00153 

2.7 

4cj1tt 

.07401 

.07362 

.00039 

.5 

IttItt 

.05453 

.05411 

.00042 

.8 

sum 

.26559 

.26199 

.00360 

1.4 

aSee  the  text  for  description  of  sets  used  and  configurations 
included. 
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2  2  * 

Computed  SO  a- 'd  Cl  expectation  values  for  the  X  and  A  i  states  of  1.10. 
All  quantities  are  given  in  a.u.  unless  explicit lv  sped  lied. 


State 

2 

2.* 

Calculation 

SCI 

Cl 

!  CF 

Ci 

Inton  uciear  distance  R 

3.20 

3.20 

5.00 

3.00 

Total  energy  E 

-e:.3U5 

-82.5342 

-82.3006 

-82.5236 

Klnatlc  energy  (T  ) 

82.296  7 

82.4999 

62 .2995 

82.5179 

Potontial  energy  <V> 

-164,6082 

-165.0341 

-164.6004 

-165.0415 

VlrUl  coefficient  <T)/<V>  -2.0002 

-2.0004 

-2.0000 

-2.0001 

<1/«u> 

8.2477 

8.2487 

6.3888 

8.3902 

<l/rQ> 

21.3975 

2  3.3898 

23.4420 

23.4361 

A 

ft* 

"Se 

1  V 

w 

v’* 

39.8157 

39 . 1107 

40.1022 

39.6198 

<ZL1> 

28.3040 

26. 2589 

26.3962 

26.  3431 

r  pole  Booent  t  (Debye) 

(LlV)  6.872 

6.758 

6.090 

5.955 

<P> 

8.4421 

8.5265 

8.8565 

8.9349 

<ru> 

31,2585 

31.2866 

29.6294 

29.6564 

<ro> 

16. 3530 

16.4630 

15.9807 

16.0926 

<ZLr. 

95.6060 

95.5425 

81.0025 

82.9119 

<P3> 

11.4507 

11.8318 

12.6056 

12.9988 

<ru> 

107.0567 

107.  3741 

95.6081 

95.9107 

Quedrupolc  oom.  Q  < 10“ 2b 

e»u  cn‘)a  5.471 

5.544 

7. 120 

7.210 

,'zu/ru> 

0.7831 

0.7789 

0.8891 

0.8655 

<  Vr0  > 

-0.2949 

-0.2914 

-0.3344 

-0.3338 

Force  on  M 

0.006 

-0.007 

0.001 

-0.010 

Forca  on  0 

-0.016 

0.011 

-O.009 

-0.004 

Not  forco  on  nuclei 

+  tT0  -0.010 

0.006 

-0.008 

-0.014 

<<3Zll  *  rU,/:rU'> 

0.7695 

0.2674 

0 . 3006 

0.2965 

<OZZQ  -  r*)/2r>. 

0.8544 

0.6589 

-1.2661 

-1.2886 

qu  -  -(.*»/«*) u 

-0  .051 

-0  .04  7 

-0.009 

-0.004 

,0  •  -<J2V/iZ2)0 

-1.526 

-1.535 

2.754 

2.799 

*0  la  computed  with  reepect  to  the  center  of  n*s 

lor  7U160. 

abl^  IV.  Molecular  properties  for  the  II  state  of  LiO  computed  with  the  Cl  (1632  configurations) 
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AE  is  difference  between  the  Cl  and  HF  energies, 
polarity  is  Li+0“. 

;  ©  is  computed  with  respect  to  the  center  of  mass  for  7Li 


2  + 

Table  V.  Molecular  properties  for  the  Z  state  of  LiO  computed  with  the  Cl  (1458  configurations) 
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Table  VI.  Computed  dissociation  energies  and  term  energies  for  LiO. 


Calculation 

HF 

Cl 

State 

2U 

Y 

2ll 

2  + 

l 

E(R  )  (a.u.) 
e 

-82.3115 

-82.3008 

-82.5342 

-82.5239 

E(R=50)  (a.u.) 

-82.2427 

-82.1756 

-82.4104 

-82.3434 

D  a  (eV) 
e 

1.87 

3.41 

3.37 

4.90 

Te(2iI-2Z)  (cm-1) 

2342 

2330 

a  2 

The  observed  value  for  D  (X  H)  is  3.39  +  0.26  eV  taken  from  Ref.  1. 
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Table  VII.  Computed  spectroscopic  constants  for  LiO.  Units  for  the 
spectroscopic  constants  are  cm 


Isotope 

State 

Calc . 

(a) 

e 

(d  X 
eAe 

v(l-0) 

B 

e 

a 

e 

y» 

7Li160 

2 

X^H 

Cl 

851.5 

12.5 

825.2 

1.202 

0.0151 

1.695 

HF 

854.4 

7.1 

840.4 

1.2  35 

0.0175 

1.673 

2  + 

AE 

CI 

866.7 

6.7 

853.4 

1.349 

0.0199 

1.599 

HF 

881.5 

8.3 

865.0 

1.376 

0.0213 

1.585 

6Li160 

x2n 

CI 

899.4 

14.2 

870.1 

1.341 

0.0178 

HF 

902.4 

7.9 

886.8 

1.378 

0.0206 

a2z+ 

CI 

915.5 

7.5 

900.5 

1.505 

0.0234 

HF 

931.0 

9.3 

912.8 

1.535 

0.0252 

/  a  yi.'ifi  i.ins 
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